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SPARGER SYSTEM FOR MMH-HELIUM VENTS 

Based on a calculated vent flow rate and MMH concentration, a Tl-59 
program was run to determine total sparger nole area for a given sparger 
inlet pressure. Hole diameter is determined from a mass transfer analysis 
in \he holding tank to achieve complete capture of MMH. In aoaition, 
based on oxidation kinetics and vapor pressure data, MMH atmospheric con- 
centrations are determined 2 ft above the holding tank. 
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A. MMH Vent Rates 

A typical venting profile is given on the attached TPS. With 800 gallons 
of ullage and the assumption of ideal bahavior (in this ''T-essure range 
helium exhibits heating rather than cooling- upon expansion, but only about 
0.06°C per atm; therefore, isothermal expansion is reasonable) the 
volumetric rate would be 1459 scfm during the first 10 seconds. With a 
vapor pressure of 0.957 psia the MMH concentration would be about 
3400 ppm. At the completion of the vent the rate would be about 162 scfm 
and the concentration (assuming rapid vapor equilibrium during venting) 
16,780 ppm. 

For the purpose of designing a sparger, let's take 1500 scfm as the max 
flow and 100 psi as the sparger inlet pressure (either an analysis of A P 
in the vent system or pressure measurement during a vent is necessary to 
determine this parameter). 

For the purpose of analysing dispersion of MMH to the atmosphere above the 
holding tank, we need to estimate the total amount of MMH vented. The 
average flow rate is 800 scfm for 3 minutes, i.e., we're venting 2400 scf 
which is 6.7 lb moles. At an average concentration of about 10,000 ppm 
MMH, we would be venting 0.067 lb moles of MMH. 


8. Sparger Program 

The attached sparger program will determine the total hole area and area 
distribution required for a liquid sparger. 

Since our case is a gas dispersed into a liquid phase, we have to employ 
the gas orifice equation instead of equation (9) in the attached article 
by Mink. 

The gas orifice equation: 

W = 1891 Y d Q 2 C * 

Where W = lb/hr 

AP = psi 
Vi= ftVlb 

d 0 = orifice diameter, in. 

converts to 

a = .26 Q . using a value of 0.6 for both Y 

1830 [ A>/p ] * 

and C. 

Therefore, equation (9) in the article by Mink probably overpredicts the 
area by a factor of 4. Keeping this in mind, the conditions of 1500 scfm 
and 100 psi were run, giving the following results: 
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For 


A = 11,220 GPM 
B = 100 psi 

C = 20 psi (includes 9-ft head in tank) 
D = 2 inches (Sparger I.D.) 

A' = 0.011 
B' = 0.011 
C' = 7 ft 


a = 0.72 in 2 


Since this is an overpredicticn, we probably only need 0.18 in 2 total 
hole area (of course, it is extremely important to know the actual sparger 
inlet pressure; our value of 100 psi is used for illustrative purposes). 


C. Mass Transfer Analysis in Holding Tank 


1. Rate 

The mass transfer rate to the tank can be approximated by: 


Where 


n a = \J 


4ir v 

^AB t 



N _ g moles transfered to water 
M CM Z sec 


D AB = diffusivity of MMH in H 2 O 
Vt = terminal velocity of rise of drop 
C a 0 = solubility of MMH in H?0 in moles/volume 
dp = drop diameter 

In order to use this equation we need to determine the drop diameter, 
dp, and the terminal velocity. 


The terminal velocity is given by: 




Where 4 is the surface tension of water and is liquid density. 
The particle (drop) diameter in inches is given by: 
d p = 0.279 (Re 0 )-°- 05 

based on a jet correlation (p. 141, Treybal). 


Using an upper limit of 50,000 on Re 0 , we get 0.16 inches for the bubble 
diameter. This gives a value of 0.77 ft/sec for Vt- 


♦ 
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Plugging these values into the equation for gives a very high rate of 
mass transfer indicating that most, if not all of the MMH, should be 
absorbed. 

2. Orifice Diameter and Number of Holes 


The orifice Reynolds's number Re 0 = 



K 


Where Vp is the velocity in the orifice and p and y are density and 
viscosity, respectively. 

Using the criterion of Re 0 = 50,000, we get the following two equations: 
For 1500 scfm 


1. n d 0 = 11.36 


2 . 



n = a = 0.18 in 2 


where n = no. of holes 

Therefore it d Q (11.36) = 0.18 
4 

or d 0 = 0.020 inches 
. * . n = 568 

So 568 holes at 0.02 inches per hole ( 1/50 ") should do the job in this 
case. 


D. Mass Transfer Analysis from Tank to Atmosphere 
This analysis is based on the following assumptions: 

1. 0.067 lb moles of MMH are completely absorbed from the sparger 
during a vent in a well mixed tank of water 10' X 5' X 9' filled to 
the 6- foot level. 

2. MMH is depleted by two mechanisms 

a. diffusion from the gas liquid interface through 3 ft of air 
to the top o ' the tank where it is swept away (zero concentration 
at top) 

b. in the tank MMH disappears by oxidation at a rate of 100 
ppm/day (This number is an average of open pan tests done by Eric 
Miller and data in the Florida Tech handbook.) 
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3. the Henry's law relationship for NMH is y = mx, where m = 0.05, 
and y and x are the mole fractions of MMH in air and H2O, 
respectively 

4. equilibrium exists at the gas liquid interface 

5. diffusion in the gas phase is quasi-steady and unidirectional 

First, let's compute the initial mole fraction of MMH in the tank, 

Xf 0 , in ppm: 

300 ft 3 H^O X 62.4 lb, X 1 lb mole = 1040 lb moles H„0 
* ft 3 18 lb 2 

.067 = 64 ppm = Xf_ 

1040 0 

The diffusion rate in the gas phase is given by 
kl «$MMH-Air A C 

A = X 


Where - flux in moles/area time 

MMH-Air - diffusivity of MMH in air 
A - cross sectional area of tank 
C - gas molar density 

X/ - mole fxn of MMH in tank at any time, t 

jP - distance between the gas-liquid interface and the top of the 
tank 

Using a value of 0.2 cm^/sec = 0.775 ft2 

hr 


foro&MMH_Air we 9 et 

N A = 0.00177 Xj* (t) moles 

Fir- 

Now we can write a mass balance for MMH in the tank to compute X^ as a 
function of time, t. 

rate of accumulation = - rate of depletion - rate at which MMH leaves 
of MMH in tank by oxidation thru interface 

or 

1040 dX £ = -10' A (1040) . 1 day - 0.00177X jt 

“dt day 24 hours 

where t is in hours 
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The solution to this equation is: 

X, = 2.35(e-l*7 x 10- 6 t_i) + X o e-1-7 X 10'h 

0 


or X^ = (Xc 0 + 2.35) (e-1-7 X 10- 6 t) 
for complete depletion X^ = 0 and 

so -1.7 X lCT 6 t 2.35 

e ' 2.35 + 64 X lCf 6 


2.35 

Xj( o = 64 X 10-6 
0.99997 


-1.7 X 10“ 6 t = -0.0000272337 
; t = 16 hours 

Lastly, let’s look at the concentration two feet above the gas-liquid 
interface initially (this will drop off to zero at 16 hours). 

The concentration profile above the liquid is: 

_ dy^ = 1.06 X 10-6 
dx 

initially at the gas liquid interface we have Y A = 0.05 (64 X 10“6) 
(where X = 0) 

Y a = 3.2 X 10-6 = 3.2 ppm 


so 3.2 - Y? ft = 1.06 
2 ft 


Y 2ft = l* 08 ppm 

Any effort to speed up the oxidation process such as low ph should enable 
the maintenance of low concentration above the tank. Forced convection 
would also do the job. 


E. Recommendations 


1. Determine the holding tank inlet pressure during a vent 

2. Measure MMH in the tank and in the vapor space above the tank 
during a vent with current system 

3. Do 1 and 2 after sparger is installed 

4. Consider agitation if bubble coalescence occurs 

5. Determine treatment schedule and/or second stage process to meet 
EPA standards for atmospheric MMH 



ORIGINAL PAGE IS 
OF POOR QUALITY 


NASA JSC WSTF 

TEST PREPARATION SHEET , 

CONTINUATION SHEET h=r=r 


DESCRIPTION CONTINUED 


ft ^ 


^ Vxi 






ft\ AO 


SjiQ. | 


W 7 


\W 1 


s 


Wo 




\\* 


\ 0 ^c 




%\ 


\ 40 | 1 % 


|Q£9WI 


!B 


ajg 

\ Sc 

— 

£S 


\Uft 

SM 


Mo 

s>i 

eh 


Mo 

S3> 

El 


\^5 




OO 


s 


'SO 


IS 


mi 


lie r msi ( r« v r«t ?i) 


19-46 




































i 


■v\ 




. I 


original page is 
OF. POOH QUALITY 


Hole-area distribution for 
liquid spargers 


This TI-59 program does the tedious calculations 
involved in designing spargers, and prints (or 
displays) the results. 



William H. Mink, BaiUlU Columbus Laboraloruj 


O In many operations, it is necessary to sparge • liquid or 
gas into another fluid; frequently, a uni form distribution is 
desired. Since, owing to frictional losses, the pressure 
decreases as the fluid flows down the sparger pipe, the hole 
area must increase along the length of the sparger to main- 
tain uniform flow from it. 

H. W. Coopert h?s developed a aeries of relationships 
that are useful in calculating sparger hole area for 
uniform distribution. It has been pointed out that these 
equations arc in error because the originator neglected 
to include the effect of velocity head. This velocity-head 
effect can be significant. Only when the pipe LfD ratio 
is high or when flowrates are low do the Cooper 
equations provide an approximation. The author has 
developed a new program, brsed on a stepwise 
calculation, which yields correct hole sizes regardless of 
L/D ratio or flowrates. 

Calculational approach 

The pressure causing flow through any orifice is the 
difference between the static pressure at the point of the 
orifice and the ambient pressure The static pressure is 
the total pressure less the velocity head and any fric- 
tional loss: 



r, * f, - A* - r, 

(<) 


For the tparger, the total inlet preuure it: 




(2) 


In Eq. (1) and (2), 


f 

_ F*P V*P 

(3) 
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The friction factor, F f depends on the Reynolds num- 
ber, 

F = I6/// iU ,ifAf iU < 2,100 (6) 

F m 0.0035 + if £ 2,100 (7) 

where 

124 rap 


AV 


(8^ 


The orifice area can be calculated by the orifice equa- 
tion, which can be simplified to: 


a = £/l830f(^ * 


(9) 


A program has been developed for the TI-59/PC- 
100 A, using the above procedure. The program calcu- 
lates the hole area for each 10% of sparger length (using 
E«. (9)). 

Ihe program is shown in Table L Certain constants 
for the calculations and for the alphanumeric opera- 
tions must be stored in the data registers. These arc 
shown in Table n. Once the calculator has been pro- 
gramed and the constants stored, both the program and 
data register comcnu may be recorded on magnetic 
cards for future use. 

User instructions for calculation of the sparger hole 
areas that will achieve uniform flow distribution are 
given in Table III. 

When running this program, the printer firm prints 
out the entered data. Neat, it prinu the hole area for 
each 10% of sparger length, and finally U prinu the 
total sparger hole area. 


Procedure 

For the calculation of the hole area in the first section 
(first 10% of length), the flowrate down the pipe is set at 
♦ •»«* *« HaI#.**** 


$ 
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Program tor uso In calculating hot* area of liquid tpargaro 
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04 

4 
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85 

4 
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93 

, 

426 

00 

0 

427 

00 

0 

428 

03 

3 

429 

05 

5 

430 

95 

■ 

431 

42 

STO 

432 

16 

16 

433 

61 

CTO 

434 

8$ 
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lation, the fluwrate down the pipe is reduced by <0% of 
♦he entered value. 

The velocity is calculated by Eq. (4) (Step 157 in the 
program shown in Table I), the velocity head by Eq. (3) 
(Step 172), and the Reynolds number by Eq. (8) (Step 
189).* The appropriate equation for friction factor is 
selected [Eq. (6) or (7)]. The frictional loss isxalculated 
by Eq. (5) (Step 222). 

The pressure difference, P t — P ft is calculated using 
Eq. (1) and the entered ambient pressure. This differ- 
ence is then used to calculate the area, using Eq. (9) 
This seiics of calculations starts at Step 251. 

It snould be noted that, since the entered P in is a 
static pressure, P ( i9 must be calculated using Eq (2). 
(Since holes are assumed to be located at the center of 
each section, the frictional lots for the first section takes 
place over only 5% of the pipe length rather than *0% 
as for the subsequent sections. To adjust for the pro- 
gram, which subtracts the frictional loss for 10% of the 
pipe length for each hole calculation, the initial tot? 1 
pressure used in the calculation is increased by one-half 
•f the initial frictional loss.) 

Sample problem 

Cooper demonstrated the calculation of hole area 
with the following example: 

Consider 2,000 gal/min of water at 25 psia and 95*F, 
which is to be sprayed onto a distributor plate via holes 
in an internal pipe that is 10 fi long. Determine the area 
allocation for a 6-in. pipe. Tower pressure is 24.5 psia 
In this cue, p xr 62.4 Ib/cu ft, p = 0.76 cP. 

Table IV shows the printout for this example. The areas 
calculated differ substantially from those calculated by 
Cooptr. in this example, velocity head it an important 
factor. 

•TV I.tmUi m **t*mm4 V I •* tMiWU O) 4mm 

mm •• mAmm |tiwm it* xi •» mtm mm 

•i im <*4 w ip **»• r ■ a iw «r 1 v» 

lAa, i im 'tor rr mli. 



Fig. I shows, diagrams tically, the application of the 
calculated results. (The hole tiac and spacing shown are 
for illustrative purposes and do not relate to the sample 
problem.) In this example, the hole area could be ob- 
tained by several combinations of hole sixes and num- 
bers of holes. One possible arrangement would be to 
have 16 holes over the first 10% of length. 1 1 l-in 
holes over the next 10% of length, and finally 6 I -in holes 
over the last 10% of length. 

If the TI-59 is to be used without the primer, it will 
be necessary to make one change in the program, ax 
follows: 

I. Enter ihe program as shown in Table 1 
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HOLE-AREA DISTRIBUTION FOP LIQUID SPARGERS 


Sample calculation Table IV 


LlOUlB Sf'hPCEF: 


2000. 

CF’H 

25. 

P-JN 

?4. 5 

P*ftH 

6. 0648 

IB 

0. 7 C 

vise 

62. 4 

BEN 

10. 

L 


2 

12. 945 

0/10 

S. (-97 

/20 

7. 101 

/$0 

6. 230 

/40 

5, 600 

/SO 

5. $07 

/60 

5.04 7 

/70 

4. see 

/$0 

4. 74 e 

-'90 

4.677 

/too 

TOTAL HOLE AREh» 

IN* 

65. 296 



2. Mike a program addition by prcoing the follow- 
ing key sequence: CTO; 344; LRN; 2nd Del; 

2nd Del; LRN; GTO; 312; LRN; 2nd Inr, R/S; 
LRN. 

With this change, the calculator will display the hole 
area for the Tint 10% of sparger length after E is pressed. 
Repeatedly pressing R/S will cause the calculator to 
display the hole area for each of the next sparger seg- 
ments in turn and, finally, the total hole area. As an aid 
to the operator, for identification purposes, the hole 
areas for the sparger segments are displayed in Fix 3 


Nomenclature 

s Orificr area, in. ? 

4 Sparger pipe 1.0,, in. 

D Sparger pipe I.D., ft 

L Pipe length, ft 

P 9 Ambient pressure (pressure sparger is discharging 
into), ptia 

P lm Initial (entered) static pressure, psia 
P ( Frictional loss, psi 

P, Static pressure, psia 

P, Total pressure, psia 

Velocity head, psi 

{£ Flowrate entering sparger pipe, gpm 
V Velocity in sparger pipe, ft/s 
p Fluid density, lb/ft* 
ft Fluid viscosity, eP 


mode while the toul hole area is displayed in Fix 9 
mode. 

Should the data entered result in an impossible solu- 
tion (such as might occur if the sparger diameter were 
too small to permit the required flow at the specified 
pressures), the calculator will print out “SELECT 
LARGER PAIN OR ID. M 


For HP-67/97 users 

The HP version of the program closely follows the Tl 
version. Table V prvides a listing of the HP program, user 
instructions for the HP version are identical to those pro- 
vided in Table III, with A\ B', and C\ replaced by a, b, 
and c, respectively. 
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on 
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31 
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-M 
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m? 

3e 05 
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KM 

nix 
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001 5 

002 5 
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6^4 fW 
025 0U5 
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OOP 
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•tor 
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•to# 
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•to# 
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Cato 
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Cato 
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2 

02 

Ml 

/X 

54 
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02 
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Oi 
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-£ 
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7 

07 

144 
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55 15 

105 

T 


-35 

10* 
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60 
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11 
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A 
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4 
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KU 
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•0 

104 

1 


OJ 
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4 

•-24 
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o; 

1*7 

J 


03 
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•35 

105 

0 


00 

120 
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run 
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H* 

1 


-35 
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35 

07 
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127 

5)00 

55 14 
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0* 
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31 
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A 


54 
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54 
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04 

no 
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2 

•i 
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Si 03 
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* 


-3: 

m 
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12 

152 

mo 
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153 

4 
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A 
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135 
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Oi 
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5* 

05 
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x* 

55 

154 

CHS 
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Oi 
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55 
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4 

04 
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;i-n 

J37 

4 


'24 

114 

r 


•35 
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CFI 
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22 OJ 

133 
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24 

n* 

7 


•: 

1J? 

3 

05 

J5P 

ttflf 
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51 
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-24 

HI 

« 

-55 

157 

0 


02 
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SPARGER SYSTEM FOR MMH-HELIUM VENTS 


On page A the Sparger program for an example vent condition was 
run giving a result cf a = 0.72 in^ (total Sparger area). It 
was decided to run the same program for the same conditions, 
except for the Sparger inlet pressure which was 100 psi 
previously and is now taken to be 21 psi (instead of a A p of 
80 psi as in the previous case we now have a4 P of 1 psi 
through the Sparger). 


The 


results are the following: 

Hole Area, in 2 / l0% of Sparger length, L %L 

. 664 10 

.605 20 

.562 ’ 30 

.529 40 

.503 50 

.484 60 

.468 70 

.457 80 

.449 90 

.444 100 


Total Area 


5.166 in 2 = a 


* 


r 
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Again, since the program is for liquids we need to multiply 
these numbers by 1/ A to get the areas needed for a gas vent . 
Also, in this case the area required continuously decreases as 
we go down the Sparger to compensate for velocity head loss. 

Using this result, if we now refer to page 5 we have the 
following two equations for the first 10% of Sparger: 

1. ndo = 1.136 (10% of flow for each 10% of 

Sparger length) 

2. TT(do) 2 n = a = . 664 = .166 in 2 

2 4 


Therefore n do (1.136) = .166 

A 


and do = .184 inches 

n = 6.17 holes 


Table I shows the results for 10% increments in length down the 
Sparger . 
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TABLE I. HOLE SIZE AND NUMBER REQUIRED 


Percent 


n (no. of holes) 


do , in 


the previous case. 


10 

6.17 

.184 

20 

6.80 

.167 

30 

7 . 2? 

.156 

40 

7.75 

.147 

50 

8.15 

. .139 

60 

8.47 

.134 

70 

8.76 

.130 

80 

8.97 

.127 

90 

K\ 
« — < 
• 

On 

.124 

100 

9.23 

.123 

we need far 

fewer holes of a larger 

diamei 


V 

t 

t 

r 
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l 
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